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^\ | Abstract. - We perform spectral simulations of dynamo for magnetic Prandtl number of one with 

Taylor-Green forcing. We observe dynamo transition through a supercritical pitchfork bifurcation. 

■_^ ■ Beyond the transition, the numerical simulations reveal complex dynamo states with windows 

(-^ , of constant, periodic, quasiperiodic, and chaotic magnetic field configurations. For some forcing 

amplitudes, multiple attractors were obtained for different initial conditions. We show that one of 



the chaotic windows follows the period-doubling route to chaos. 



The phenomenon of spontaneous generation of magnetic fields in magnetohydrodynamics 

\Q ' (MHD) is known as the dynamo effect. This is believed to be the generating mechanism for 

the magnetic field in astrophysical bodies such as planets and stars pQ. Dynamo has been 

observed in several numerical simulations and laboratory experiments in which a magnetic 

field is amplified significantly from a seed magnetic field. 

Some of the important problems in dynamo are related to the dynamo transition, e.g., 
what is the nature of the bifurcation from the fluid state to the dynamo state with a magnetic 
field?, how is chaos generated in dynamo?, what leads to the reversals of magnetic field?, how 
is dynamo affected by forcing and geometry?, etc. In this letter we attempt to understand 
the nature of bifurcations and the routes to chaos in dynamo under the Taylor-Green (TG) 
K> , forcing. 

!-h ' The governing equations for dynamo are 



d t v + (v-V)v = -Vp+(B- V)B + iA7 2 v + F, (1) 

<9 t B + (v-V)B = (B-V)v + r/V 2 B, (2) 

V-v = 0, (3) 

V-B = 0, (4) 

where v is the velocity, B is the magnetic field, p is the total pressure (thermal+magnetic), F 
is the external force field, v is the kinematic viscosity, and r\ is the magnetic diffusivity. Some 
of the important parameters that characterise dynamo are the magnetic Prandtl number 
Pm = v/r], the Reynolds number Re = UL/v, and the magnetic Reynolds number Rm = 
UL/rj, where U and L are the characteristic velocity and length of the system respectively. 
Only two of the above three parameters are independent (Rm = Re * Pm). Note that 
the range of magnetic Prandtl number Pm observed in nature is extremely large. Liquid 
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metals and solar plasma typically have small Pm (10~ 5 to 10 -2 ), while interstellar medium 
has typically large Pm (of the order of 10 14 ). There are other important parameters in 
dynamo, for example, Rossby number, Ekman number, Raylcigh number, etc. To abstract 
and simplify the physics of dynamo transition, in the present letter we focus only on basic 
MHD equations [Eqs. (QUI)] that have only two independent parameters, Pm and Re. 

In MHD systems, forcing is typically applied only to the velocity field. A steady mag- 
netic field is generated beyond a certain critical forcing. The magnetic Reynolds number 
corresponding to this forcing is called the critical magnetic Reynolds number {Rm). The 
first bifurcation is generally a pitchfork bifurcation. Subsequent secondary bifurcations yield 
more complex configurations like periodic, quasiperiodic, and chaotic magnetic fields. In this 
letter we will investigate the bifurcations and routes to chaos for Pm = 1 under the TG 
forcing. 

Dynamo transition has been studied extensively using experiments, numerical simula- 
tions, and theoretical modelling. The Riga dynamo experiment [3J, which relies on the 
concept of the Ponomarcnko dynamo [3], was the first successful liquid metal experiment 
in which dynamo was observed. Dynamo has also been observed in the Karlsruhe exper- 
iment [4] in which the flow pattern is helical. Recently Monchaux et al. [5j performed 
dynamo experiments on liquid sodium (Pm ~ 10 ) confined within a horizontal cylinder. 
The fluid was forced using two fans at the two ends of the cylinder. The speeds of the 
two fans were used as the main controlling parameters. This experiment is called VKS (von 
Karman Sodium) due to the nature of its forcing and the velocity configurations. Monchaux 
et al. [5] observed dynamo for Rm ~ 30 for equal and opposite forcing frequencies of the 
two fans as a result of a supercritical pitchfork bifurcation. The resulting magnetic field 
is axisymmctric and is constant in time. By changing the relative frequencies of the fans, 
the VKS experimental team could obtain a variety of dynamo solutions including periodic, 
quasiperiodic, and chaotic magnetic field configurations [BJ. They also observed reversals of 
the magnetic field. 

Various researchers have simulated dynamo using direct numerical simulations (DNS) 
with random, TG, ABC, and Roberts forcing and observed dynamo in both low and high 
Pm regimes. Variations of R C M as a function of Pm have been investigated in these studies. 
TG forcing has certain similarity with the von Karman flow configuration and the VKS 
experiment, hence it has become quite popular. Nore et al. [7] demonstrated the existence 
of dynamo action under TG forcing. Ponty et al. [8j applied TG forcing and reported that 
R C M increases sharply with P M X as turbulence sets in, and then saturates. Ponty et al. [9] 
observed subcritical dynamo transition by changing Rm in their simulations. Mininni |10j 
observed various dynamical regimes including time-periodic oscillations and well-defined 
spatial structures. Dubrullc et al. [llj reported various bifurcations in TG flows for both 
hydrodynamics and MHD simulations. 

In another set of simulations, Schckochihin et al. [12| and Iskakov et al. [13] applied non- 
helical random forcing and showed the existence of dynamo for both large Pm and small Pm ■ 
Podvigina [14] . Mininni [T5], and Mininni and Montgomery [TB] simulated helical dynamo 
using ABC and Roberts forcing. Podvigina P3] simulated dynamo with ABC forcing and 
studied various magnetic field states including chaos. The author related these states with 
inherent symmetries of the system. Gissinger et al. [17j generated dynamo by axisymmetric 
forced flow in a spherical domain; they also observed chaotic magnetic field reversals in the 
same geometry. Morin and Dormy [18] studied dynamo in a rotating spherical shell and 
observed either supercritical or subcritical dynamo transition depending on the chosen set 
of parameter values. Glatzmier and Roberts [TO] simulated geodynamo and observed several 
interesting phenomena including the reversals of the magnetic field. 

Dynamo transitions have also been studied using low-dimensional models that are con- 
structed through heuristic arguments, or from partial differential equations using Galerkin 
projections. These models are amenable to analytical investigations of bifurcations due to 
smaller number of modes. Rikitake model or disk dynamo, which belongs to the former 
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category, is a discrete model with conducting discs and current-carrying wires [20] . This 
model shows complex dynamical behaviour including constant, periodic, quasiperiodic, and 
chaotic magnetic fields. The model of Petrelis and coworkers [H] is based on amplitude 
equations and symmetry arguments. They attempted to understand the origin of various 
dynamo states and magnetic field reversals using the saddle-node bifurcation in their low- 
dimensional model 21 . Models based on Galerkin projections have been constructed by 
Donner et al. [22] and Verma et al. [23] . The magnetic field in the six-mode model of Vcrma 
et al. |23| is generated by a supercritical pitchfork bifurcation at the transition point. This 
model exhibits only a constant magnetic field (in time), and it does not show complex dy- 
namical behaviour. On the contrary, the 152-mode model of Donner et al. [52] shows various 
patterns including constant, time-periodic, quasiperiodic, and chaotic states for Pm = 1. 
There are other dynamo models based on scale separations, e.g., a-dynamo, /3-dynamo, 
etc. [J. 

In this letter we numerically study the dynamo transition for TG forcing given by 



F(fco) = F 



sin(kox) cos(fcoy) cos(fcoz) 

— cos(kox) sm(koy) cos(fcoz) 





(5) 



where Fq is the forcing amplitude, and fco is the wavenumber of the forcing taken as 2 in 
this work. Note that the helicity of the force, (VxF).F = everywhere. Nore et al. [7] and 
Ponty et al. [5] however argue that local fluctuations in kinetic helicity are generated by the 
above forcing. We solve Eqs. (HE]) using TARANG [H], a pseudo-spectral code, in a 3D 
box of dimensions 2n on each side for Pm = 1- We apply Runge-Kutta fourth order scheme 
for time advancement. The time increment dt is determined using the CFL condition (dt = 
Ax/(20y/E u ), where Ax is the grid size, and E u is the total kinetic energy). The number of 
grid points used in our simulation is 64 3 . Our runs are dealiased using 2/3 rule. The range 
of Reynolds number investivated is from 6 to 160, for which our simulations are well resolved 
as k m axV (the largest wavenumber times the Kolmogorov length) is always greater than 1.3. 
This observation is corroborated from the well resolved kinetic energy spectrum (E u (k)) 
and magnetic energy spectrum (E b (k)) shown in fig. Q] for v = 77 = 0.1, and F = 4.8. The 
number of interacting Fourier modes of our dynamo system is 64 3 , which is quite large. Our 
simulations however reveal that only a small fraction of them carry most of the energy. In 
many numerical runs with different Fq's, we observe that the most prominent velocity Fourier 
modes are (±2, ±2, ±2), (±4, ±4, ±4), (±4, ±4,0), (0,±8,±4), and the most prominent 
magnetic Fourier modes are (0, 0, ±1), (0, 0, ±2), (0, 0, ±3), (±2, ±2, =p3), (=p2, ±2, ±1). Here 
the three arguments refer to x, y, and z components of the wavenumber. Note that other 
modes like (±4, 0, ±4) for velocity field are also present due to symmetry. The most energetic 
velocity Fourier mode is (±2, ±2, ±2) due to fco = 2 of the TG forcing. Among the magnetic 
modes, the most dominant modes are B(0,0, 1) and B(0, 0,2). The mode B(0,0, 1) is 
generated due to the nonlinear interactions between (v(2, 2,2),B(— 2, — 2, — 1)), and the 
mode B(0, 0,2) is generated by (v(2, 2, 2), B(— 2, — 2, 0)). We observe a dynamic interplay 
between B(0,0, 1) and B(0, 0, 2) modes. The other important magnetic mode is B(0,0,3) 
which participates with (v(2, 2, 2), B(-2, -2, -1)). 

We perform DNS for various forcing parameter Fq. Initial transients are discarded and 
only the steady-state configurations are analysed. For F = 4.8, the steady-state velocity and 
magnetic fields (snapshots of the magnitudes) are shown in figs. [5Ja,b) respectively. These 
figures indicate that the TG forcing yields well defined velocity and magnetic structures. 
For fco = 2, chosen for our runs, the simulation box has 16 TG cells for the velocity field. 
As shown in fig. [2^b), the magnetic energy is concentrated in two major slabs along with 
two minor slabs. We also observe that \B Z \ <C \B X \, \B y \. The above configuration is due to 
the prominence of the modes B(0, 0, 1) and B(0, 0, 2). Note that the z components for these 
modes are zero due to the incompressibility condition. Thus, we can understand the global 
TG structures in terms of the dominant Fourier modes. Figure [3] contains the vector plot 
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of the magnetic field (x and y components only) for a cross-section of fig. [2jb) at z = 4.08. 
The magnetic field in this plane is approximately along —45°. These field configurations are 
similar to those reported by Mininni et al. |f0j . 

To explore various bifurcations of the dynamo state, we vary Fq in the interval [f : 40]. 
We observe pure fluid solutions (E h = 0) till Fq = 3.9, after which nonzero steady-state 
magnetic field (E b > 0) emerges with each component of the magnetic field as constant (in 
time). The magnetic Reynolds number at this transition regime is approximately 19.9, hence 
R C M ~ 19.9. The number of variables of our dynamo system is unfortunately rather large 
(~ 64 3 ). We could however focus on the most energetic modes that determine the system 
dynamics. We investigate the time series of these modes, and obtain various dynamo states: 
constant (fixed point), periodic, quasiperiodic, and chaotic magnetic fields. In fig.E]we show 
the time series of a fixed point (F = 3.9), a periodic state (F = 10), a quasiperiodic state 
(F = 36), and a chaotic state (F = 4.8) of the dynamo. 

Different states of a dynamical system are elegantly illustrated in the bifurcation diagram 
that contains information about the birth of new states. In fig. [3] we construct a bifurca- 
tion diagram by plotting the time averaged value of the amplitude of the magnetic mode 
(5(0,0,1)) for different Fq (obtained from around 60 DNS runs). Nonzero magnetic field 
appears at around Fq = 3.9 through a pitchfork bifurcation. There is no hysteresis near the 
onset of dynamo as Fq is increased or decreased. Hence, the dynamo transition is through a 
supercritical pitchfork bifurcation. After the primary instability (or bifurcation) , we observe 
different kinds of dynamo states like constant, periodic, quasiperiodic, and chaotic magnetic 
fields as evident from the bifurcation diagram. The windows of these states appear for var- 
ious range of F , e.g., chaos appears for F = 11 — 12. We also find windows of F where 
(5(0,0,1)) becomes relatively small, but (5(0,0,2)) becomes significant. The dynamics of 
interchange of energy between B(0, 0, 1) and B(0, 0, 2) is not apparent at present, and it will 
be studied in future. 

A careful analysis of the dynamics for a given forcing reveals coexistence of multiple 
dynamo states. We illustrate this feature using state space or phase space. The phase 
space of our dynamo is very large (~ 64 3 ). However, projections of the phase space on 
the subspace formed by the most energetic modes contain most of the information of the 
system. For Fq = 16, two different sets of initial conditions yield either a fixed point or a 
periodic solution. These two states are illustrated in figs. EJa,b) using projections of the 
phase space on (5(0, 0, l)-5(0, 0, 3)) plane. At F = 36, the two different coexisting dynamo 
states involve quasiperiodic and chaotic magnetic fields as illustrated in figs.[6](c,d). 

The bifurcation diagram (fig. [5]) is quite complex, and it is not possible to probe all 
the secondary bifurcations using DNS. In the present letter we focus on a narrow window 
(Fq = 4.6 — 4.8) of fig. [S]to investigate how chaos appears in this window. As we will describe 
below, here we observe a period-doubling route to chaos. To understand the transition to 
chaos better, we plot the phase space projections on the (5(0,0, l)-5(0, 0, 3)) plane. At 
around Fq = 4.6 we observe a fixed point or a constant magnetic field as evident from 
fig-EKa). At Fq ~ 4.73, the fixed point bifurcates to a periodic solution (limit cycle) through 
a Hopf bifurcation (also shown in fig. [3(a)). The power spectral density (PSD) plot of 
the time series for this Fq exhibits a peak at a single frequency /i = / ~ 0.018 (in non- 
dimensional units) as shown in the right side of fig. Eta). As we increase Fq, we observe 
period-2, period-4, and period-8 solutions at F = 4.77,4.777, and 4.778 respectively. These 
new states are generated through the "period-doubling bifurcations" . Projections of the 
phase space on the (5(0, 0, l)-5(0, 0, 3)) plane along with the power spectra for these three 
states are shown in figs. [3Jb,c,d) respectively. The power spectra clearly indicate peaks at 
subharmonics fi = //2, f 3 — //4, / 4 = //8 corresponding to the period-2, period-4, and 
period-8 dynamo states. At Fq ~ 4.8, the system becomes chaotic as evident from the 
phase space projection and the broad frequency spectrum shown in fig. [7£e). The above 
observations indicate that chaos appears at Fq ~ 4.8 through the period-doubling route to 
chaos. The value of magnetic Reynolds number for these states is around 22. The bifurcation 
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diagram (fig. [5]) exhibits several other windows of chaos whose origin has not been explored 
in this letter. 

In conclusion, our numerical simulations of dynamo for Pm = 1 with Taylor-Green 
forcing reveal that the dynamo transition takes place through a supercritical pitchfork bi- 
furcation. After the primary bifurcation, the system exhibits several windows of constant, 
periodic, quasiperiodic and chaotic solutions. Wc also find multiple coexisting attractors for 
a given parameter; different initial conditions take the system to one or the other attractor. 
A careful analysis of one of the chaotic windows reveals that the dynamo becomes chaotic 
through a period-doubling route to chaos. 

Our numerical simulations reveal several dynamo states that have been observed in 
experiments (such as VKS), earlier numerical simulations, and low-dimensional models. The 
geometry and the forcing of our simulations however are simpler than those of experiments. 
Yet the above similarities may be due to certain inherent common features of dynamo. 
Future numerical simulations with more realistic geometry and forcing functions will reveal 
valuable insights into this puzzle. 

We thank S. Fauvc, E. Dormy, D. Carati, K. Kumar, and T. Lcssinnes for fruitful 
discussions and comments. This work was supported by a research grant of DST India as 
Swarnajayanti fellowship to MKV. 
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Figure captions 

Fig. [1] The kinetic energy spectrum (E u {k)) and the magnetic energy spectrum (E h (k)) for 
Fa = 4.8 for which the magnetic field is chaotic. The spectra indicates that the simulations 
are well resolved. 

Fig. [2J Flow structure of the Taylor-Green (TG) flow: (a) volume rendering of the velocity 
field magnitude that illustrates 16 TG cells due to fco = 2; (b) volume rendering of the 
magnetic field magnitude illustrating two major and two minor slabs. This emphasizes the 
dominance of B(0, 0, 1) and B(0, 0, 2) magnetic modes. 

Fig. [3j The magnetic field vector plot for a cross-section of fig. EJb) at z = 4.08. 

Fig. |U The time series of the real part of B x (0, 0, 1) for F = 3.9(a), 10(b), 36(c) and 4.8(d). 

Fig. [5j Bifurcation diagram exhibiting various dynamo states of our simulations with FP = 
fixed point, P = periodic state, QP = quasipcriodic state, and C = chaotic state. The time- 
averaged amplitude of 5(0, 0, 1) is plotted for various Fq. Multiple attractors are observed 
at Fo = 16 and Fq = 36 (see fig.E]). The circled region of the inset exhibits period-doubling 
route to chaos (see fig. [7]). 

Fig. [5J Phase space projections on the (-B(0, 0, 1)-_B(0, 0, 3)) plane for Fo = 16(a,b) and Fo = 
36(c,d) exhibiting multiple attractors. Figures (a,b) show coexisting fixed point and periodic 
states for two different initial conditions. Figures (c,d) show coexisting quasipcriodic and 
chaotic states. 

Fig. [71 Phase space projections on the (B(0, 0, l)-5(0,0, 3)) plane and the power spectra 
exhibiting period-doubling route to chaos: (a) a fixed point at Fo = 4.6, and pcriod-1 
solution at Fo = 4.73; the power spectrum in RHS exhibits a peak at / ~ 0.018 (non- 
dimensional units) for Fo = 4.73; (b) period-2 solution at Fo = 4.77 with power spectrum 
showing two peaks at fa and fa{= //2); (c) period-4 solution at Fo = 4.777 with power 
spectrum with three peaks at fa, fa and fa(= //4); (d) period-8 solution at Fo = 4.778 with 
power spectrum showing four peaks at fa, fa, fa and fi{= //8); (e) chaotic solution with a 
broad-band power spectrum. 
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